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One of the most important challenges for the largely accepted idea that Galactic CRs are acceler¬ 
ated in SNR shocks is the maximum energy at which particles can be accelerated. The resonant 
streaming instability, long invoked for magnetic field amplification at shocks, can not provide suf¬ 
ficiently high fields and efficient enough scattering so as to ensure particle acceleration up to the 
knee. Here we discuss the non-resonant version of this instability which, with its faster growth and 
larger value of the amplified field, increases the achievable maximum energy. Because of their 
higher explosion rate, we focus on type II SNe expanding in their red supergiant wind and we 
find that the transition between Ejecta Dominated (ED) and Sedov-Taylor (ST) phases takes place 
at very early times. In this environment, the accelerated particle spectrum shows no high energy 
exponential cut-off but a spectral break at the maximum energy (Em). Moreover, the maximum 
energy of protons can easily reach PeV energies. With this model, we tried to fit KASCADE- 
Grande and ARGO -YBJ data but failed to find a parameter combination that can explain both 
data sets. We discuss the different scenarios implied by the two data sets. 
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1. Introduction 


Supernova remnants (SNRs) are considered the most probable sources of Galactic CRs based 
on energetic argument s, Supporting evidence has been found in the X-ray band where amplified 
magnetic fields have been inferred at the shock with surrounding medium (for reviews see, e.g., 
10 0 )- Magnetic field amplification, not only downstream but also upstream of the shock, is 
fundamental in order to obtain high energies. The quest for understanding what kind of instability 
could lead to the right conditions is still open. The resonant instability induced by CR streaming 
leads to estimate Em ~ 10^ — 10"^ GeV still well below the knee\ this is due to its resonant 

nature that fixes the saturation at 8B/B 
studied in Refs. 


1 122, 38, ra, 33, 341. Non-resonant modes, largely 


and [10], starting with a wavenumber k much larger than the particle gyro- 
radius but then forming fluctuations on larger scales allowing resonant scattering, grow in a very 
short time. This kind of instability may allow to reach the knee if applied to the case of supemovae 
expanding in the wind of their red super-giant (RSG) progenitor star [0,0], where the highest 
energies are reached within a few decades after the SN explosion, namely before the beginning of 
the ST phase of the explosion. 

The origin of the knee of the all-particle spectrum is a very challenging and a very important is¬ 
sue because it is strongly correlated with the transition between galactic and extragalactic CRs. In 
Ref. ]^, in order to fit the Auger data S0I.0. that require a predominant light component at 
EeV energies, the authors introduce an extra component of protons; this added proton population 
is in good agreement with the proton spectrum measured by KASCADE-Grande [||] (see also the 
results of ICETOP [§]). However, proton and iron spectra detected by KASCADE-Grande in the 
energy region lO'^ — lO'^ eV make the issue of understanding what the maximum achievable en¬ 
ergy is even more challenging: these data can be explained only with a CR spectrum that is not 
cut-off exponentially at the maximum energy, or with a new unknown, class of sources with maxi¬ 
mum energy much in excess of the knee. Eurthermore, recently, ARGO-YBJ ]|T^ and YACl-Tibet 
Array ]^] experiments measured a knee in the light component at ~ 650 TeV, appreciably below 
the PeV energy of the all-particle knee, making the situation ever more complex. 

In this work we analyse the probability of reaching the energy of the ’knee’ in particle accel¬ 
eration during the expansion of a SN shock in the dense region occupied by the wind of a RSG. In 
particular, we investigated in detail the implications of the so called Non-Resonant Hybrid (NRH) 
instability described by Refs. ] [T^ , by computing the maximum energy and the overall 

particle spectrum produced during the whole SN expansion, both in the ED and ST phases, and 


comparing it with KASCADE-Grande and ARGO data ]]15[]. 


2. Model 


2.1 Maximum Energy 

The non-resonant instability has a very high growth rate and, despite the fastest mode has a 
wavelength much shorter than the Earmor radius of the particles generating the current, an ’inverse 


cascade’ mechanism during its non-linear evolution, increases the spatial scale involved [10, 29, 
14]. Particles with energy equal to the current achievable maximum, Em, move from downstream to 
upstream where, by definition, they can not be scattered resonantly and escape the system moving 
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quasi-balistically at a speed close to that of light. In order to obtain an estimation of the maximum 
energy achievable, we passed through different steps. First of all, we write the current, at distance 
R from the explosion center, due to particles escaped when the shock location was at r: 

jcR{R,r) = ncK.x{R^RM{r))eVsh{r) = (^) (^. 1 ) 

where ^CR is the CR acceleration efficiency and 




{Em/Eq) In {Em/Eq) 


P \Eo) 


1 -te) 
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i8/0 


( 2 . 2 ) 


The function accounts for normalisation of the particle distribution function which is taken 

to be fs{E) oc and extending between a minimum energy Eq, that does not depend on time 

and a maximum energy Em which does depend on time, as we will see (details in [|T^]). Numerical 
simulations show that the saturation of the instability occurs after N, ~ 5 e- folds [|T^], namely 


jMdt = J {kMVA)dt 


■ 5. 


(2.3) 


where va is the Alfven speed in the unperturbed magnetic field Bq. The wavenumber of fhe fasfesl 
growing modes, kM, can be easily esfimafed considering balance befween currenf and magnefic 
fension: 

471 

— — jcR, (2-4) 

_ c 

We can see fhaf if depends on CR currenf (Eq. |2.1[ ) and, consequenfly, on fhe densify profile. For 
a supernova expanding in fhe uniform ISM, fhe densify can be assumed consfanf p{R) = PisM- 
For a fype II supernova expanding ifs pre-SN wind, instead, fhe densify profile can be written as 
p{R) = ^ ~ P{Ro) (x)^’ where M is fhe rate of mass loss of fhe red gianf and is 

fhe wind velocify. In fhe following we will simply wrife p{R) R^’” and consider fhe fwo cases: 
m = 0, fhaf corresponds fo a uniform medium, and m = 2, fhaf describes expansion in fhe progenifor 


wind. Af fhis poinf, differenfiafing Eq. 2.3 wifh respecf fo R, we obfain an implicif expression for 
fhe maximum energy: 

2e 
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-^CRVs{Rf^4np{R)R^ . 


(2.5) 
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From fhis equafion, we can hnd maximum energy for fype I and fype II SNae for differenl injecfion 
specfral indices. In fhe case of a E^^ (jS = 0) source specfrum, for fype I evenfs we can esfimafe: 
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It is apparent that for standard parameters pertaining the two types of explosions, type II SNe can 
reach about one order of magnitude larger maximum energies than type la, and in particular that 
if particles are accelerated and escape as described here, these sources easily seem to be able to 
provide proton acceleration up to the knee (details in [[Isl]) . 


2.2 Particle spectrum 

In order to obtain the particle spectral behavior provided by NRH instability, we use the ex¬ 


pression in Eq. 2.1 for the current and assuming a general power-law dependence of the shock 
radius on time: R oc which also implies oc So, the escaping particle spectrum is: 


N,. 


AE) = 

e dE 


4?rgcR 

EoS{X,m] 


pv]R^X{E) , 


(2.8) 


where x{E) = {d/dE)[\ /*P(E)) is a function of E alone, to be computed from Eq. and S{X,m) 
depends on the details of the SNR evolution (more detailed explanation can be found in [jlsl]). 
When computing x{E) explicitly from Eq. we find: 
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where the sign ps refers to the fact that we have used the assumption E » Eq. It is apparent that 
the power-law dependence of x{E) can only be —2 or steeper. Therefore the spectrum of CRs 
released during the Sedov-Taylor phase is the same as the spectrum of accelerated particles in the 
source only if the latter is E^'^ or steeper, while for flatter source spectra Nesc{E) E^^. 

In the following, we take the value of A appropriate to describe the different evolutionary 
stages of the SNR from Ref. [ jl^ [T^ , ^n[ These works give the remnant expansion law during 
the different stages for generic density profiles of both the SN ejecta oc R ^ and the ambient 
medium p oc 7?^'”. Once A is given, from Eq. ^ and ^ we can derive the spectrum released by 
the SN during the ED and ST phases as [|I5|]: 


Nesc{E) < 


Nesc{E) < 


E ED phase; 

m = 0,k = 7 

£■-( 2 + 6 ) 

£-(4+3e) EQ 

m = 2,k = 9 

£-(2+e) 3E pj^^gg. 


( 2 . 10 ) 


( 2 . 11 ) 


This result is very interesting because it shows that the spectrum of escaping CRs integrated over 
time during the SN shock expansion is a broken power law, with an index close to 2 at low energies 
(belowEM) and steeper at higher energies. There is no trace of the exponential cut-off expected 
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by standard DSA theory. The transition between the two power laws at the maximum energy Em 
corresponds to the transition between the end of the ED phase, where most mass has already been 
processed, and the beginning of the ST (adiabatic) phase. 

In the next step we take into account the transport of nuclei in the Galaxy. Assuming that 
the transport has reached a stationary regime and that it is dominated by diffusion and that both 
injection and spallation occur in a thin disc of size 2h, the solution of the transport equation in the 
Galactic disc can be easily written as follows: 


hAp) 


Nescjp)^ X{p) 
kRI^v 1 + ^’ 


( 2 . 12 ) 


where we introduced the grammage A(p) = Ud^vnip-^^. Here is the gas density of the Galactic 
disc, assumed of thickness 2h and radius Rd, and the escaping spectrum Nesc{p) is calculated as 
discussed above but using a parametric form for Rsh in order to model the transition from the ED to 
the ST phase []T5|]. Written in this form, the solution is very clear. At particle momenta for which 
the grammage is small compared with Xa the standard solution is recovered fo,aip) = 
while in the range of momenta at which the solution is spallation-dominated, the observed spectrum 
reproduces the shape of the injection spectrum. 

In order to calculate the spectra of CRs at the Earth, an evaluation of the grammage is needed 
and we use the formulation of Ref. [^| as a starting point. Given that the CR spectrum observed 
at Earth has an energy dependence derived from TRACER and CREAM data [^, ^ ], 

in our calculations we used a slope of the diffusion coefficient 5 = 2.65 — pinj, where for we 
considered the values 2 and 2 . 31 . In the first case, the diffusion coefficient has a slightly stronger 
energy dependence than the one used by |^] A in the other case, instead, we use the exact form 
given in that work. Eor the spallation cross section we adopt the simple formulation of Ref. 12C ] 
(details in [p3l]). 


3. Results 


Type II SNe are the most likely sources to accelerate protons up to the knee. We focus on this 
type of explosions in order to hnd a ht for the all particle spectrum. With this aim we vary the 
different parameters, considering the right ejecta density prohle. We have seen that the value of 
Em is determined by the combination of the CR acceleration efficiency, i§cr> and the energetics of 
the SN, Esn- On the other hand, the flux of CRs at the Earth derives from a combination of 
Esn and the rate of SNe, 9^. Eitting the proton spectrum in the knee region imposes an additional 
constraint on the combination Esn, (see Eig. |^. 

We calculated the required efficiency and rates necessary to reach the knee and ht the overall 
CR hux: for steep spectra at the shock, the escape current is lower and this forces one to have 
larger SN energetics and larger CR acceleration efficiencies, which is counterintuitive for a steep 
spectrum of accelerated particles. Eor hat spectra at the shock (EA, the spectrum of escaping 
particles is also hat and PeV energies can be reached for ordinary parameters of the SN explosion. 
Reaching PeV energies is even easier for hard spectra at the shock. Our calculations show that a 
type II SN with standard energetics (Esn = 10^' erg, ^CR = 10%) can accelerate CRs up to energy 
of the knee, Em ~ 1 PeV. The spectrum obtained with these parameters is shown in Pig. 3 of 115]. 


5 









NRH instability and type II Supemovae 


Martina Cardillo 



Figure 1: Relations between physical parameters in the case of injection index p = —2. The shaded area 
in the ^CR — Esn plane indicates the allowed range of parameters to reach Em = 1 PeV. The different curves 
indicate the combination of values of the parameters for which the observed proton flux is also fitted. Each 
line refers to a given SN rate (as labeled). 

Expressing the time dependence of the maximum energy through the parametrization used in 
this work, we can show that PeV energies can be reached at very early time, ~ 10 — 30 years after 
the SN explosion, . This result not only changes the SNR paradigm according to which the highest 
energies are reached several hundred years after explosion, but also narrows down the probability 
of catching a PeVatron in action in our own Galaxy. 

Finally, we used our model for an attempt to fit recent challenging data of KASCADE-Grande 
and ARGO. We can fit individually both sets of data, using reasonable values of the SN parameters 
and CR acceleration efficiencies; by fhe way, fhere is no model fhaf can fil bofh af fhe same lime (see 
Fig.|). KASCADE-Grande dala require a SN energy Esn = 2 x 10^'erg, an efficiency ^CR = 20% 
and a rale of explosion of = 1/110 yr^' lhal lead lo a maximum energy Em = 3.7 PeV, whereas 
we can fil ARGO dala wilh Esn = 10^' erg, 91 = 1/15 and ^CR — 5.2% lhal lead lo Em = 507 
TeV 




Figure 2: (Left) our best fit model for KASCADE-Grande data [|^], with k = 9, Esn = 2 x 10^' erg, 91 = 
1/110 yr^', Ej^ = 3.7 PeV and ^cR — 20%. The highest energy data are fitted with an "ad hoc" extragalactic 
component. (Right) Best fit model to ARGO data [ [Tsl [l^ , with k = 9, Esn = 10^* erg, 91 = 1 /15 yr^*, 
Em = 507 TeV and ^CR = 5.2%. 


The different results obtained by these two experiments measuring the same quantities points 


6 












NRH instability and type II Supernovae 


Martina Cardillo 


to the possible underestimate of systematie uneertainties involved in the adopted experimental teeh- 

niques. 
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